Multi-particle correlations of an oscillating scatterer 
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Using multi-particle distribution functions we calculate the correlations produced by a periodi- 
cally driven scatterer in a system of noninteracting electrons at zero temperature. The multi-particle 
correlations due to a quantum exchange effect are expressed in terms of photon-assisted scattering 
amplitudes. The results we obtain are valid for slow but arbitrary in strength driving. We show 
that even for large amplitude pumps the zero-frequency noise power is related to two-particle corre- 
lations. In addition to two-particle correlations a large amplitude pump can generate multi-particle 
correlations. 
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I. INTRODUCTION 

In recent papersii- it was shown that a periodically 
driven unbiased scatterer, a quantum pump^'^'^'^, can gen- 
erate entanglement in a system of noninteracting fermions. 
Entanglemenl!^*^*— *— is a property of quantum objects which 
has a powerful potential for applications in information 
processing,^^ The main resource for quantum computation 
are specific correlations between different parts of a quantum 
system. These correlations reflect the fact that some quan- 
tum states of a whole system can not be represented as a 
product of states corresponding to its different parts. Such 
states are called non-separable. For them the joint probabil- 
ity to find the system's parts in some states does not factorize 
into the product of probabilities to find any of the parts in 
a corresponding state. Thus the joint probability consists of 
two terms: the first one is a product of probabilities for single 
parts, and the second one depends on correlations between 
the system's parts. The last term in itself can be used as an 
entanglement measure (see, e.g., Refe.fS 14}. Such a measure 
is especially appropriate for systems of non-interacting par- 
ticles where there are no other sources of correlations except 
the quantum-statistical interactions due to the Pauli exclu- 
sion principle. 

In open (i.e., coupled to external reservoirs) solid state sys- 
tems it is more appropriate to discuss the properties of states 
rather then the properties of separate individual particles. 
The occupation of a single-particle state is given by the single 
particle distribution function which is a probability to find 
the particle in a state under consideration. By analogy, for a 
system of noninteracting particles, one can consider the prob- 
ability to find several particles in some given (single-particle) 
states. This probability is a multi-particle distribution func- 
tion, which is a diagonal element of a corresponding multi- 
particle density matrix. If the particles are correlated (i.e., 
if the multi-particle state is not-separable) then the multi- 
particle distribution function does not factorize into a product 
of single-particle distribution functions. 

In this paper we calculate the multi-particle distribution 
functions for particle flows generated by a periodically driven 
scatterer and relate them to the noise power which can be 



measured in solid state structures. We consider a mesoscopic 
scatterer coupled to Nr stationary reservoirs of noninteracting 
electrons via single channel leads. We will number electron 
states in leads with Greek letters, a, (5, ... , which include both 
the orbital (a lead number) and the spin indices. We suppose 
that all the reservoirs are at zero temperature and have the 
same Fermi energy ^. The electrons in different reservoirs are 
uncorrelated. Therefore, the multi-particle incoming state is 
pure and separable. After scattering by the pump the state 
remains pure.- However it becomes correlated and thus non- 
separable r-»2i Correlations between the out-going particles ap- 
pear as a result of an interplay of two factors, the Pauli ex- 
clusion principle and photon-assisted scattering at a working 
pump. 

Below, we consider, flrst, the effect of photon-assisted scat- 
tering which appears already in single-particle scattering. 
Then, in a second step, we analyze two- and multi-particle 
correlations. 

The paper is organized as follows. In Sec|n]we discuss the 
approximations made and analyze the single-particle distribu- 
tion function for electrons scattered by the pump. In the next 
section, Sec lIIII we calculate the two-particle correlations and 
their relation to the zero-frequency noise power in the case of 
large amplitude pumps. Then, in Sec liVI we consider multi- 
particle correlations produced by the pump. We conclude in 
Sec El 



II. SINGLE-PARTICLE SCATTERING 

In the presence of a scatterer oscillating with frequency S7, 
particles absorb or emit energy quanta while traversing the 
pump. In contrast, while in the leads, they are described by 
the stationary Schrodinger equation. Thus one can classify 
the states of electrons (both out-going and incoming) in the 
leads according to a particle energy E like in the stationary 
case. Therefore, we will use the second quantization operators 
corresponding to states of particles with definite energy. We 
denote the annihilation operator for incoming states aa{E) 
and the one for scattered electrons be (E). These operators are 
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related to each other through the Floquet scattering matrix 
Sf of a periodically driven scatterer^^*— *— 



A. Adiabatic approximation 

If the pump oscillates slowly, S7 — > 0, then the Floquet scat- 
tering matrix can be expressed in terms of the Fourier coef- 
ficients Sm{E) of the stationary scattering matrix S{E, {pi}) 
with time dependent parameters pi{t + T) — Pi{t), where 
T = 2Tv/n, as foUowspi^ 



SF{E,Err.)^S-„,iE) + 0{n). 



(1) 



Here Em — E + mhQ; m is an integer. We ignore terms of 
order Q, and higher. The Floquet scattering matrix couples 
electron states at energies shifted by one or several quanta 
hO,. In what follows, we will consider only electrons with 
energy close to the Fermi energy /i. Thus it is convenient to 
introduce the following notation for energy variables: E = 
tn = e + nhil, where the Floquet energy e lies within the 
interval, fi < e < fi + h^l. As a consequence, the operators for 
out-going particles relate to the ones for incoming particles 
through: 



1 (tn) — Sa 



fH jn — m 



ap{€„ 



(2) 



We evaluate the stationary scattering matrix at the Fermi en- 
ergy. Current conservation forces the Floquet scattering ma- 
trix to be unitary. In the adiabatic approximation of Eq.Q 
the unitarity reads: 



lOOafi, 



(3) 

This follows directly from the fact that the stationary scat- 
tering matrix is unitary: S = SS^ = T. Here X is a unit 
matrix. 

At zero temperature all incoming states with em < fJ, 
(m < — 1) are filled while the ones with em > fJ- (m > 0) are 
empty. The corresponding single-particle distribution func- 
tion fa"\em) = {ciL{em)aa{em)) IS the Fermi distribution 
function which is the Heaviside step function at zero temper- 
ature. 



B. Non-equilibrium out-going particles 

The situation is different for out-going particles. Due to 
interaction with an oscillating scatterer an electron can gain 
or lose some energy quanta h^l. Therefore, some out-going 
states with e„ > can be partially occupied while the ones 
with en < II become partially emptied. The probability to 
find an out-going particle in the state with energy e„ in lead 
a is given by the single-particle distribution functioni— 

— oo 

/i°""(en) = (6l(6„)6.(fn))=^ |S„^,„-™|'. (4) 

7 m — — 1 

Here {• • ■ ) denotes a quantum-statistical average. Since the 
6-operators are expressed in terms of a-operators [see, Eq.l|5J] 



we, in fact, average over (the product of) the equilibrium 
states of reservoirs. 

To assess at energies far from the Fermi energy we 

take into account the following. For each particular pump, 
there exists Umax such that S±n ~ for jn| > Umax- Then far 
below (above) the Fermi energy, n < —Umax {n > +nmax), all 
the relevant incoming states are filled (empty) . Therefore, the 
out-going state is filled (empty) with a unit probability. For 
n < —Umax this follows directly from the unitarity condition 
Eq.@ taken at / = and a — /3. 

In contrast, in the vicinity of the Fermi energy, e„ « /i (with 
1^1 < nmax), the distribution function Eq.(0J describes par- 
tially occupied states above and below fi (see Fig0. Since 
for these energies not all the relevant [i.e., the states with 
energy em for which -Umax < n — m < +nmax] incoming 
states are filled (empty). The incoming states correspond- 
ing to n — Umax < JTi < — 1 are filled while those with 
< m < n + Umax are empty. Close to the Fermi en- 
ergy particles in UmaxNr incoming channels are scattered into 
"irimaxNr available out-going channels. Since the number of 
available out-going channels is larger then the number of in- 
coming channels, the appearance of a particle in some out- 
going state is a random process. This results in a partial oc- 
cupation of the corresponding single-particle states and leads 
to the appearance of particle-particle correlations. 

Notice, the pump is able to correlate only those out-going 
particles which have energies shifted by one or several quanta 
fiQ,. While the particles having different Floquet energies, 
e 7^ e', remain uncorrelated like the incoming particles. 

To investigate the correlations we calculate the joint prob- 
ability to find several out-going channels occupied and com- 
pare it to the product of occupation probabilities of individ- 
ual channels. The single-channel occupation probability is 
the one-particle distribution function. In analogy we calcu- 
late multichannel occupation probabilities which are multi- 
particle distribution functions. First we consider the two- 
particle probability. 
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FIG. 1: Occupation of single particle states e„ = e -f nhO. 
in incoming channels (a) and in out-going channels (/3) of a 
working pump at zero temperature. The filled box area is 
proportional to the occupation probability. The case shown 



corresponds to n„ 



= 3, i.e. S, 



I3a,±\n\ 



0, for n > 4. /i is 



the Fermi energy; is a driving frequency. 
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III. TWO-PARTICLE SCATTERING 

To calculate the joint probability to register one out-going 
particle with energy e„ at lead a and other one with energy 
Em at lead /3 we introduce the two-particle distribution func- 
tion /^°^*''(en, Em) for out-going particles. This distribution 
function is the quantum-statistical average 

= (-Bj„,E„(en,e,7i)-BQ,0(en,e77i)>, (5) 
of the two-particle operator 

Ba,l3{tn,em) — bc,{en)bl3{e,n). (6) 

Using Eqs.ll^J and Q the distribution function can be ex- 
pressed as a sum of two contributions, 

one which factorizes and an irreducible part, a two-particle 
correlation function, 

(8) 

1 P=-l 

where ifa^(e„,em) = (pfji^rajbai^n)) is a matrix element of 
the pair correlatori- 

) is non- vanishing the two-particle probability 
does not factorize into the product of single-particle ones and 
we conclude that the particles are correlated. 

The electron-electron correlation function 5f^^'^\en,6m) 
characterizes correlations between the out-going particles. 
Since Sfj^^^^^ is the (negative of a) square of some quantity 
we can interpret it as a probability of some scattering process 
which is responsible for two-particle correlations. This pro- 
cess consists of the creation and scattering of electron-hole 
pairs. 



A. Two-particle correlations: an electron-hole pair 

view 

The matrix element Kapi^n, em) of a pair correlator is a 
sum of particular amplitudes 

^a,f3 i^ri^^m) ~ Saj.n — pSp^^m—p^ (9) 

which can be related to the excitation of an electron-hole 
pair from the filled state with energy ep incident from reser- 
voir 7. Formally one can consider the hole (an empty state) 
as a particle with corresponding second quantized operators 
a^^^ = a\ = 6^. Therefore, scattering of holes is described 
by the adjoint scattering matrix 5^. We define holes for any 
energies, above and below the Fermi energy. The hole can be 
scattered by the pump like an electron. Both electron and 
hole can emit/absorb some energy quanta hfl. Therefore, the 
amplitude Eq.@ corresponds to the process with a final state 
in which an electron with energy nhQ is scattered into lead 
a and a hole with energy mh^l is scattered into lead 13. To 
get the whole amplitude ifa/3(en,em) it is necessary to sum 
over all the indistinguishable processes leading to a given final 
state. Different processes are possible due to the presence of 



many leads (index 7) from which the pair is born and due to 
different filled energy channels (index p). Alternatively the 

quantity (^A''^'^\e„, em)^ is an amplitude for electron and 
hole to be scattered to leads /3 and a, respectively. Therefore, 
the two-particle correlation function & }'^^*'\en, tm), Eq.l|HJ, is 
the probability to create an electron-hole pair and to emit its 
constituents (an electron and a hole) into two chosen states. 
Notice, the shot noise generated by the pump is due to the 

19 

same processes^ 

We remark that the electron-hole pair shot noise generated 
by a microwave field applied to the contacts of the sample 
has recently been measured. '^^ In this experiment there is no 
app lied dc-voltage and hence no dc-current. Refs. |^ and 
I22I propose experiments with multiple microwave fields of the 
same frequency but with possible phase-lags applied to dif- 
ferent contacts of the sample. In contrast to the discussion 
provided here, since the microwave fields are applied to the 
contacts, the scattering matrix is (apart from self-consistent 
effects^^) stationary. 

One can ask why the correlations between two electrons de- 
pend on processes involving electron-hole pairs. That is due 
to the fact that in each state the (operators for the) numbers 
of electrons and holes sum up to unity. Therefore, the fiuctua- 
tion of the number of electrons and holes are not independent 
from each other. As a result the various two-particle correla- 
tion functions involving electrons and/or holes are related to 
each other in a simple way: 

'5/^7'(^n,e„) = -5/l7>(0.„,e„) 

(10) 

= -5/^T'(^n,0.„J=5/(7'(0.„,0.„). 

Here Ot„ means that there is no particle (there is a hole) in 
the state e„. Note, we define an electron-hole distribution 
function, say /^°^''(en, Ot^), by Eq.© with an electron op- 
erator 6/3 (em) being replaced by the corresponding hole oper- 
ator 6^''' (e m) = bl{em). The various two-particle distribution 
functions sum up to unity 

Ao^t)( f \ I f(o"*)cn f \ 

Ja.0 \'^ri,tm) ^ J ^0 l,OE„,tm; 

(11) 

This normalization condition justifies the point of view that 
is a two-particle joint probability. In addition Ea. lllH 
shows that the electron-hole processes affect the electron- 
electron probability. 

Note, that the whole electron-hole scattering amplitude 
Kc,i3{€n,em) is a sum of a number of two-particle (electron- 
hole) photon-assisted scattering amplitudes ^[7/3 ' (^"' 
which interfere between themselves. When all the relevant 
incoming states are filled then owing to current conserva- 
tion [i.e., to a unitarity of scattering, Eq.lO] the electron- hole 
amplitudes cancel each other and the correlations disappear. 
This holds for particles with energies far below the Fermi en- 
ergy. That is quite similar to the stationary case when the 
unitary scattering process does not produce additional corre- 
lations if incoming states are not correlatedf^^*— In contrast 
for the particles with energies close to the Fermi energy not 
all the relevant incoming state are filled and the correlations 
appear. 
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B. Pauli exclusion principle and two-particle joint 
probability 

The correlations under consideration are quantum- 
mechanical because they originate from the Pauli exclusion 
principle. According to this principle two particles can not 
be simultaneously scattered into the same out-going state, 
and any incoming state can not be a source for two parti- 
cles. Formally the Pauli principle leads to the existence of 
two amplitudes (direct and exchange) describing scattering 
of two fermions. In our case this becomes evident if we ex- 
press the two-particle distribution function in terms of Slater 
determinants: 



-oo — oc 



7 5 



-lq = -l > 



(12) 



^y^(2) / 0(x^ ,n — p SaS,n — q 



Thus, the joint probability to measure two particles in some 
given out-going states la; e„) and em) is a sum of various 
probabilities describing photon-assisted scattering of two elec- 
trons. Each such scattering is described by the Slater determi- 
nant whose elements are photon-assisted single-particle scat- 
tering amplitudes. The sum of squared direct/exchange two- 
particle amplitudes gives the product of two single-particle 
distribution functions /a°"*'(en)/^°"'^(em) (it is a joint prob- 
ability for uncorrelated particles). While the interference 
between direct and exchange amplitudes is responsible for 
electron-electron correlations. 



C. Current noise and two-particle correlations 

The quantity 5/„°^''(en, tm) characterizes how much two 
particles in the states |q; tn) and £m) are correlated be- 
tween themselves. To measure these correlations (as well as 
to measure single-particle distribution functions) one needs to 
use energy filters with resolution better then hfl. Apparently 
this is a challenge for an experimental realization especially 
in the adiabatic regime, f2 — > 0. Therefore, it is useful to 
relate what we calculated to what can be easier measured ex- 
perimentally. In mesoscopic systems the zero frequency noise 
power is an observable which characterizes the correlations 
between the particlesi2S*2Si2i 

The zero frequency noise power, 

T oo 

— OO 

_ (13) 

produced by the pump was calculated in Ref.0. Here ia{t) is 
the current operator in lead a and AIa{t) = Ia{t) — {ia{t)). In 
terms of (5/^°^'' (e„, e™) the zero-temperature current-current 
cross- correlator (a 7^ /3) reads 



the statistically independent sets of states corresponding to 
the different Floquet energies < e — ^ < HQ. 

We see that the zero-frequency noise power depends on 
pair correlations between the out-going particles in different 
energy channels exiting the pump through two chosen leads 
a and 13. Therefore, it can be used as a measure of two- 
particle correlations produced by the pump. That agrees with 
conclusions made in Ref. However the pump produces the 
multi-particle correlations as well. 



IV. MULTI-PARTICLE SCATTERING 



To show this we introduce an A'^-particle operator 



J = l 



and the corresponding A'^-particle distribution function 

fiT-,aK = (Sqi, •■ ,c«jv)i which Can be represented 

in terms oi N x N Slater determinants 



r ) M<. E 



N'. ^ E 

71...7JV P1...PN = -1 

I,--- ,iV. 



r(out) 
/cK 1 , ■ ■ ■ 

(16) 

The A^-particle correlations can be defined as follows. Let 
us subtract from the A'^-particle distribution function var- 
ious products of distribution functions and correlations 
functions corresponding to a smaller number of parti- 
cles. Then the remaining A^-particle correlation function 



by A'' scattered particles. In terms of elements of the pair 
correlator it reads as follows: 



will describe the correlations shared 



sfLi,-- ,Qjv (^"1 1 ■ ■ ■ 5 ) — (1) E 

(ir) 



Here the sum runs over the set of all non- equivalent permu- 
tations P^r ~ (ri,r2, • • ■ ,rjv) of integers from 1 to N. Non- 
equivalent permutations of integers are those that can not be 
obtained via a cyclic permutation of each other. 

Notice, the multi-particle correlation function, Eg. 1)17^ . for 
N > 2 can be either positive or negative unlike the two- 
particle correlation function, Eq.JSJ, which is definitely nega- 
tive. 



a/3 



(out) / 



0, c^l3. (14) 



n— — oc m— - 



A. Three-particle correlations 



The auto-correlator P^o, can be calculated using the conser- 
vation law "Ylia ~ 0. Notice the factor Q/2n counts all 



As an example we consider the three-particle distribution 
function in detail (for brevity we suppressed corresponding 
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B. Generating function 
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FIG. 2: Three-particle correlation function 

_J as a function of the phase 
lag ip. The scattering potential is: V{x,t) — 
5{x + a /2) [Vo + 2Vi cos{Qt)] + S{x~a/2) [Vo + 2Vi cos{Qt + (p)] 
The parameters are: a = lOOvr; Vo ~ 20; Vi = 10; and the 
Fermi energy fi — 1.0197. We use the units 2me — h — e = 1, 
where me is the mass of an electron. 



Formally one can represent the N-particle correlation func- 
tion, Eg. 1171 . as the N-th order derivative of some generat- 
ing function dependent on the pair correlator K. To this 
end we introduce a diagonal matrix A — diag(Ap) of aux- 
ihary fields Ap., where p = 1,2, •■■ ,Nmax with Nmax = 
{2nmax + l)Nr. The index p = (a,n) includes both the 
lead index a — 1, 2, • • ■ ,Nr and the energy channel index 
n = 0, ±1, • ■ ■ , ±n,nax- Then the N-particle correlation func- 
tion can be calculated as 



crr(out) _ 
"J PI,--- ,Pn 



d^Trlnll + AK 



9Apj ■ ■ ■ dXp^ 



(19) 



Here 7 is a unit matrix of dimension Nmax x Nmax ■ For A'^ = 1 
we obtain the single-particle distribution function, Eq.(|ljl. 

To prove that Ea. H19|l leads to Ea. H17|l it is necessary to 
use the Taylor expansion of Eg. 1)19^ in powers of A. Only the 
N-th term of this expansion does contribute to the quantity 
of interest. 

We emphasize that Eg. 1191 determines a multi-particle cor- 
relation function which is a different quantity then transferred 
charge which is usually discussed;"^ The important difference 
is that here we are interesting in energy-resolved correlation 
characteristics of pumped particles rather then in statistical 
properties of the whole current carried by the particles at all 
the energies. 



energy arguments): 



/(out) 



r(out) r(out) r{out) . r{out) c r{out) 



Ja jg Jl 



r{out) <■ fiout) r (out) r /.(out) , f; f(out) 
+ .1/3 OJari + Jl "J a, 13 +"Ja,l3,j' 



C. Higher order current cumulants 

The multi-particle correlations generated by the pump can 
be related to higher order current cumulants in full analogy 
with how the two-particle correlations are related to the zero- 
frequency noise power, Eg. ll4ll . 

We consider the zero-frequency Fourier transform of the 
Nth-order current correlation function symmetrized in lead 
indices 



^fafM^"" ' ' "^"T ) = 2-Re 



X r=-l 



-rS' 



5 p=-l 



V 9=-! 



(18) 

where the single-particle distribution function /a""'' and 
the two-particle correlation function Sfl^^^^ are defined 
in Eq.Q and Eq.®, respectively. In the station- 
ary case a three-particle correlation function is positive, 

^/aX-T^*^"" ' '^"|9' "^"t) ^ 25c«,/3<5/3,7^n<,,n^5n^,n.^. lu COUtraSt 

in the driven case the sign can be either positive or negative. 



In Fig|5| we give the three-particle correlation function 
(5/{^°^'^(enQ, , £„_ J ) for a two-terminal ID scatterer com- 
posed by two oscillating 5-function potentials placed a dis- 
tance a from one another. The indices L and R stand for the 
left and right terminals, respectively. The variation of the 
phase lag ip between oscillating potentials changes the value 
and the sign of the three-particle correlation function. 



1 

aT 



-Pjv «=i 

(20) 

Here the sum runs over the set of all the permutations Pn = 
(ri, ■ ■ ■ , r-]v) of integers from 1 to N. We suppose ri = 0. 

At zero temperature the current cross-correlator (oi 7^ 
• • ■ 7^ ajv) can be expressed in terms of the N-particle cor- 
relation functions for outgoing particles as follows: 



2lT 



(21) 

Such a relation follows straightforwardly from two observa- 
tions. First, at zero temperature only the currents carried by 
the outgoing particles contribute to the current correlation 
function. And, second, the operators for outgoing particles in 
the expression for the current cross-correlation function can 
be arranged into the same order as in the expression for the 
N-particle correlation function. 



6 



Therefore, the N-particle correlations generated by the 
pump can be experimentaUy probed via the Nth-order cross- 
correlator of currents flowing into the leads attached to the 
pump. 



V. DISCUSSION AND CONCLUSION 

We have shown that an adiabatic quantum pump gener- 
ates multi-particle correlations which originate from the Pauli 
exclusion principle and appear in photon-assisted scattering 
channels. Formally the pump produces multi-particle corre- 
lations up to an infinite order (A^ oo). This is because the 
pump can excite all the incoming particles. In fact, as we 
mentioned already, this number can be bounded by rimaxNr. 

The multi-particle correlations of all orders are important 
to characterize the whole outgoing state produced by the 
pump. For instance, one projects the whole state onto a 
state with an exact number of excited particles. Such a pro- 
jected state is a multi-particle state, therefore, its statistical- 
correlation properties depend on multi-particle correlations 
as well. In Ref. Q the projected electron-hole pair state pro- 
duced by the pump was considered and the electron-hole en- 
tanglement entropy was calculated. Since the entanglement 
entropy depends on multi-particle correlations, it, in a gen- 
eral case, can not be related to charge noise representing only 
two-particle correlations, Ea. (ll4ll . However in some particular 
cases the pump produces effectively only two-particle correla- 



tions. This is the case for a weak amplitude {rimax ~ 1), spin- 
independent pump with two single-channel leads (A''^ = 2). 
In this case entanglement entropy can be related to charge 

1 2 

noise. ' ■ 

The multi-particle correlations Eg. 1171 . produced by the 
pump depend on the elements of the Floquet scattering ma- 
trix. These elements are amplitudes of scattering with ab- 
sorption or emission of one or several modulation quanta 
hQ. The maximum number Umax of energy quanta which can 
be absorbed/emitted during the scattering process define the 
maximum order of multi-particle correlations produced by the 
pump. The number rimax depends strongly on the strength 
of driving. At weak driving rimax = 1 while at strong driving 
rimax ^ 1. On the other hand the value of a multi-particle 
correlation function can be changed by changing the parame- 
ters of driving, for instance, the phase difference between two 
driving parameters. Therefore, the order and the amount of 
correlations generated by a pump can be simply manipulated 
by changing the parameters of a drive. 
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